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We develop a new method of computing radiation reaction for a point particle interacting with 
massless scalar and vector fields in curved space-time. It is based on the analysis of field Green's 
functions with both points lying on the particle world-line and does not require integration of the 
field stresses outside the world line as was used in the DeWitt-Brehme approach, thus leading to a 
substantial simplification of the problem. We start with space-time of an arbitrary dimension and 
show that the Hadamard expansion of the massless scalar and vector Green's functions contain only 
integer inverse powers of the Synge world function in even dimensions and only half- integer in the odd 
dimensions. The even-dimensional case then is treated in detail. We analyze divergencies, calculate 
higher-derivative counterterms, and find a recurrent formula for the local parts of the reaction 
force in neighboring dimensions. Higher-dimensional curved space counterterms are not simply the 
covariant generalizations of the flat ones, but contain additional curvature-dependent terms. We 
illustrate our formalism in four and six dimensions. In the first case we rederive the results of 
DeWitt-Brehme-Hobbs in a simpler way, in the second case we give a covariant generalization of 
the Kosyakov equation. The local part of the reaction force is found to contain a term proportional 
to the Riemann tensor which is absent in four dimensions. 
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I. INTRODUCTION 

Though already Abraham [l[ gave a correct formula 
for radiation reaction force in classical electrodynam- 
ics 

mi a = ^e 2 (z a ~z 2 z a ), (1) 

this problem, especially its generalization for curved 
space-time and gravitational radiation still attract at- 
tention. Physical status of the Eq. (P} also remains 
subject of discussion. In 1938 Dirac gave a consistent 
derivation of the Eq. ([1]) for a radiating charged which 
is now commonly referred as the Lorentz-Dirac (LD) 
equation Q. There exist a wide literature devoted 
to different aspects of this equation, its solutions and 
applications (see pi-[Tlj]. and refs. therein). Physical 
interpretation of the LD equation is obscured by the 
presence of the third derivative Schott term 

/Schott — Z ■ (2) 

This term can be found [l3| as improvement of the 
radiation recoil term equal to the loss of the particle 
momentum due to radiation, 

/emit — — g e ^ Z , (3) 

to ensure the orthogonality of the reaction force to 
the particle four-velocity. But, somewhat unexpect- 
edly, the Schott term leads to disturbing modification 
of the particle momentum balance and creates the self- 
acceleration problem, for which reason the validity of 
the LD equation is questioned [131 ] . A widespread 
opinion is that the resolution of these problems re- 
quires quantum theory. 
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Meanwhile, correct physical interpretation of the 
Schott term solves this controversy within the clas- 
sical theory. As was emphasized long ago by Teit- 
elboim Q, and more recently reconsidered (with all 
calculations explicit) in our paper the Schott 

term must be treated as a part of the electromag- 
netic momentum, but not as a part of the particle 
momentum. The LD equation (and its generalization 
which we will discuss here) describes a composite ob- 
ject whose energy-momentum contains the mechan- 
ical part and the Coulomb field part bound to the 
charge, which varies once particle is accelerated. The 
energy-momentum balance can be verified globally for 
world-lines which asymptotically are purely geodesic 
(non-accelerated). The Coulomb coat varies with time 
when the charge experiences an acceleration, more- 
over, as a result of interaction with the charge, the 
coat energy can be transformed to that of radiation. 
This provides a natural explanation of the Born para- 
dox. 

The LD equation was covariantly generalized to 
curved space-time with an arbitrary metric in 1960 
by DeWitt and Brehme [ijj and completed with the 
Ricci term by Hobbs [HI (for a recent review see [Hj]). 
In this case, in addition to the local Abraham force 
|T|) , there is a tail integral over the past history of the 
particle: 

mz- = eF« t0 z? + \e 2 (z a -z 2 i«) + t (jfy* 

T 

+R lP z-<zPz a )+e 2 zP J f a 01 (r,T')P(r')dr', (4) 

— oo 

where F" x t0 is the external electromagnetic field and 
f a j 3 1 ( T , T ') is some two-point tensor. It could be 
expected that the curved space LD equation which 
takes into account non-local effects, should violate the 
equivalence principle. But in the DeWitt-Brehme- 
Hobbs (DWBH) equation this violation is minimal: 
only the Ricci term is present in the local part of the 
LD force, not the Riemann tensor. In empty space 
the local part of the DWBH equation is just the co- 
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variantization of the fiat-space LD equation. The tail 
term of course makes the reaction problem essentially 
non-local, but this term is due to the global properties 
of space-time. 

Similar equation was also derived for the scalar ra- 
diation reaction fl~7l Il8j: 

+R^^z a )-q 2 f G>™ t [z (r) , z(t')) dr' , (5) 

J — oo 

where G rc t is the retarded solution of the wave equa- 
tion and m(r) = mo — q<f)(z) is the variable effective 
mass. 

The derivation of these generalizations of the LD 
equation to curved space-time in the existing liter- 
ature is based on lengthy calculations of the inte- 
gral contribution of the bound and radiative energy- 
momentum fluxes within the world-tube surrounding 
the particle world line. Meanwhile, the resulting equa- 
tion (including the tail term) depends only on the 
quantities localized on the world-line. This indicates 
that a simpler derivation should be possible which 
does not require considering the field outside the world 
line (in flat space such a calculation is known for a long 
time [3J ) . Recently we have shown that similar deriva- 
tion (restricted to four space-time dimensions) can be 
performed in curved space-time as well [l9[ . Note that 
we discuss here only radiation of non-gravitational na- 
ture, for gravitational waves the situation is likely to 
be more complicated. 

The purpose of the present paper is to generalize the 
approach of [19] to higher-dimensional space-times. 
Generalization to higher dimensions in the flat space- 
time was discussed recently in (20l - [22l | in view of an in- 
terest to cosmological models with large extra dimen- 
sions. Here we consider the curved space of dimension 
D ^ 4. Note that the cases of odd and even dimen- 
sions are essentially different: in odd dimensions the 
tail terms appear already in the flat space due to viola- 
tion of the Huygens' principle. (This effect can be eas- 
ily interpreted considering D-dimcnsional point parti- 
cles as parallel strings in D + 1 dimensions [H|.) We 
will present some general considerations valid in any 
dimensions and then specialize to an even-dimensional 
case. The mostly plus space-time signature is used. 



II. GREEN FUNCTIONS 

As in D — 4 the curved space Green's func- 
tions for massless fields in arbitrary dimensions can 
be constructed starting with the Hadamard solution. 



A. Hadamard expansion 

We start with the scalar Hadamard 0, [24| Green's 
function Gr{x, x') which is a solution to the homoge- 
neous scalar wave equation 

a x G H (x,x') = 0, (6) 



where □ = g^V^'V 1 ' is the curved space scalar 
D'Alembert operator. The procedure consists in ex- 
panding Gh(x,x') in terms of the Synge world func- 
tion 

si 

<y{x,x') = i(sx - s ) J g al3 z a z^dT, (7) 

so 

where an integral is taken along the geodesic x = 
z(t), z(tq) — x, z{t\) = x' connecting points x and x' 
and it is assumed that there is a unique such geodesic 
(which makes the construction essentially local). The 
sign is chosen such that a be negative for time-like 
geodesies. The gradients of the Synge two-point func- 
tion 

a^ = da{x,x')/dx", a a = da(x, x')/dx' a (8) 

satisfy the Hamilton- Jacobi equation 

g^Wv = g a ^ a a fj = 2a, (9) 

where (using DeWitt-Brchme conventions) we denote 
by the initial Greek letters a,/3, ... tensor indices 
transforming at x' , and by the older Greek letters 
fi, v, . . . the indices transforming at x. In the coin- 
cidence limit x — > x' (denoted by square brackets) 
one has 

[a] = lim a(x,x')=0, [a a ] = 0, [a J = 0. (10) 

A mixed covariant derivative V Q V M a(a;, x'), where ac- 
cording to the convention acts on x and V a oni', 
defines the van Vleck determinant 

*(~~')- dCt[VaV ^ (x ' X ' )] (11) 
y/g{x)g(x') 

Higher covariant derivatives of the Synge function in 
the coincidence limit are expressed via the metric ten- 
sor and the Riemann tensor, this is used to find Tay- 
lor's expansions of the two-point tensors [14] . 

Another useful two-point tensor is the bivector of 
parallel transport, which can be expressed through 
the vielbein associated with metric: g tiy {x) = 

g, a (x,x>)=e™(x)e a a (x>)6 a m . (12) 

This operator transforms indices of the type a into 
indices of the type /i, in particular, 

flp"^, x')a a (x, x') = -ap(x, x'), (13) 

where it was taken into account that the two gradients 
have opposite directions. 

The world function has a dimension of the length 
squared. In flat space-time 

a{x,x')= l -(x~x'Y. (14) 

A two-point bitensor (in particular, a biscalar) can be 
expanded in powers of Synge's world function and its 
gradients. It is worth noting that this type of expan- 
sion is different from the Taylor's covariant expansion 
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of a two-point bitensor. In particular, vanishing of 
cr(x,x') does not necessarily mean x = x' . 

Let us write an appropriate expansion for the 
Hadamard's two-point scalar function Gh{x,x') 
which is singular in the coincidence limit. For D = 4 
the Hadamard expansion contains two terms singu- 
lar in a, namely, <7 _1 and lner, the remaining part 
being regular in er. In higher dimensions one has to 
add other singular terms, and by dimensionality it is 
easy to guess that each dimension introduces an addi- 
tional factor ct^ 1 / 2 . Thus, the Hadamard expansion in 
D = 2d dimensions (d ^ 3/2 is integer or half- integer) 
generically must read 



Gh{x, x 1 ) 



1 



(27r) d 



' D 

.71=1 



^cr 1 "™ 72 +V lncr 



, (15) 



where g n = g n (x,x'), v = v(x,x') are coefficient two- 
point scalar functions. The logarithmic term, as will 
be clear shortly, is present only in even dimensions. 
Let us show, that in odd dimensions we actually have 
only odd powers of cr -1 / 2 , and in even dimensions — 
only even powers, that is, an expansion in terms of 
inverse integer powers of a. 

Substituting (fT5|) into (|6|), in the leading singular 
order we will have: 



2a^d fl g D +g D (aa- D) = 0. 
Using the relation 



Ua = D 



A 



we arrive at the equation 

a^2d^g D ~g D d, l lnA)=0, 
which is equivalent to 



(16) 



(17) 



(18) 



2d fi g D -g D d ll mA = g^ - 



2a 



K, (19) 



where h v {x) is an arbitrary regular vector field. The 
left hand side of this equation is analytic in the co- 
incidence limit, while the right hand side is not, un- 
less h„{x) = 0, so with the normalization condition 
[go = 1] we obtain 

g D = A 1 / 2 . (20) 

In the next singular order we obtain: 

2d^g D ^ 1 + go-xUa - (D - l)g D -i = 0. (21) 

For D = 3, this equation holds not for the go -it but 
for v. The Eq. (|2"Tj) has a solution 



BD-i - C ^TJI> 



.A 1 / 2 
V 1 



(22) 



which does not satisfy the required analyticity for g n 
unless C — 0, and hence 



9d-i = 0. 



(23) 



For D — ■■ 3 this means the absence of the logarith- 
mic term. Similarly, considering the equation for 
9(D-i-2k), k E N we find 



9D-i-2k = 0. 



(24) 



This means that for an even-dimensional space-time 
the Hadamard Green's function contains only integer 
negative powers of a plus logarithm and a regular 
part, while in the odd-dimensional case - only half- 
integer powers of a plus a regular part. 



B. 



Dimensional recurrent relation 



For the sequence of Green's functions in the flat 
space-time, the one in D + 2 dimension is proportional 
to the derivative of the Green's function in the twice 
preceding dimension D. In particular, in even dimen- 
sions the symmetric Green's function is the derivative 
of the order d — 2 of the delta-function: 

G D ~6 d - 2 (-a), a = (x-x') 2 /2 (25) 

and thus, G D+2 ~ dG D jda. Applying a convenient 
regularization of the delta function, 

8{{x-x'f) = Km 6(\(x-x') 2 \~e 2 ), (26) 

e— >+0 



we obtain 



G D+2 ~ dG u jde A . 



(27) 



This relation has a consequence that the Laurent ex- 
pansion of the Lorentz-Dirac force in terms of e in 
the even-dimensional Minkowski space has only odd 
negative powers, and no even terms. 

Now consider the Hadamard function in the curved 
space. Taking into account the Eq. (|2"4"1) . we change 
notation for g n in (|15[) in even dimensions as follows: 



G 



/d-2 
1 / u k 



H 



(2-k)° 



■v In a + w 



(28) 



\k=0 



where uq = A 1 / 2 and we also denote v = u^-x, w = 
u c i- Applying the curved space operator □ with re- 
spect to x M , we obtain from OGjj — the following 
equation: 



d-2 ( Du 



d — 3 / Ou\ 



T d-2 

1 



- 2(9^! - iti^ lnwo)^ -2ttij 
2{d ll u 2 + u 2 d f _ l lnu )cr M - 4u 2 ) + 



(a Ud -3 - 2d IL u d - 2 (T> J - + (4 - Ua)u d - 2 



1 



+ - {Uud-2 + 2{d^u d -i - 2u d ^id^mu )(j* 1 + 
a 

+{D - 2)u d -x + o-dUd) + lno-Dud-i = 0. (29) 

Each term in the right hand side must vanish indepen- 
dently, which gives the system of recurrent differential 
equations for Ui(x, x') . Integrating them along the 
geodesic connecting two points x, x 1 , one can uniquely 
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extract u±(x,x') through uq(x,x'). Furthermore, 1*2 
is determined through u\, etc. 

Now we want to compare uf > (x,x') in different di- 
mensions D. First of all we observe that uq = A 1 ' 2 
in any dimension. It is worth noting, however, that in 
the expansion in terms of a 



jg = A 1/2 = 1 + 1/12 R a0 a a a 



(30) 



the tensor indices a, j3 will run different range of val- 
ues. Keeping this subtlety in mind, we can write 



(31) 



for any D, D' . Now, equating the first line in (TSUI) 
to zero, we see that the product u\{d — 2) is subject 
to the same master equation in any dimension, so for 
D,D' >■ 5 one has 



D 



d-2 



- -tf_2- (32) 
The second line in (|2"9l ) gives similarly for D,D' ^ 7: 

u d' = u »^-\ (33) 
2 2 d< - 2 d< - 3 ' 

If D' = D + 2, this relation can be rewritten as 
n ,D+2 _ n d — 3 



1 



For u® +2 we obtain 



D+2 



d-l 



(34) 



(35) 



Comparing the equation for the a~ l term in D 
dimensions and the equations for the a~ 2 term in 
D' = D + 2 dimensions 



,D+2 
l d-l 



l d-l~ 



1 



1 



(36) 



Denote the part of the Hadamard function (|28| with- 
out the logarithmic term as "direct" (it will corre- 
spond to the non-scattered propagation of waves) 



d-2 



r - 1 Uk 

Hdil ~ (27T) d ^ rrd-l-k- 



(37) 



fc=0 



Differentiating Gh with respect to a with the rule 
that the functions Uk,v,w are not differentiated over 
a we find 



(27T) 



,8G 



H 



da 



E-^(rf-i-*) + - = 

fc=0 

E^s^- 1 -*)- ( 38 ) 

fc=0 



Taking into account the relation (|36|) . one obtains 



-(2*)-«7g = £5ti£ a 



d - 1 da 



k 



D 

d'-2 



fc=0 

d'-2 D > 



E 

k=0 



(2TrfGf D \ (39) 



T d'-k-l 



where D' = D + 2. 



C. Retarded Green's functions 

To define the retarded and advanced Green's func- 
tion one has to specify boundary conditions. For 
a general space-time this is not suggestive, so usu- 
ally one uses a quasilocal definition proposed by De- 
Witt and Brehme [l4j in the case of four dimensions. 
We follow the same approach here for arbitrary D. 
The prescription consists in constructing the Feyn- 
man's propagator Gf by means of a shift of a in the 
Hadamard expansion onto the upper complex plane 
a — > a + iO and then separating real and imaginary 
parts 



Gf = G(i) — 2iG so if . 



(40) 



This prescription works both for odd and even D, here 
we will concentrate on D even. In this case the rela- 
tions 



(a + ioy 



V- 



1 



(-1) 



n-1 • 



a" (n — 1)! 
\n(a + iO) = In |<j| + m0(-o) 



(41) 
(42) 



are useful and we find for the symmetric (self) Green's 
function: 



Gseif — — — \toGf — T77. — r 
2 4(2tt) 



d-l 



E 

,m=0 



{-l) m u d ^ m 8^ m \a) 



v9(-a) 



(43) 



To construct the retarded and advanced Green's func- 
tions one applies the formal global definition of the 
Heaviside function 6[E(x), x'] which is equal to one if 
x' is in the past of any space-like hypersurface £(#) 
containing x, and zero otherwise. Then the symmetric 
Green's function is split into the sum of the retarded 
and advanced ones like in the flat space: 

G rct (x,x') = 29{E(x),x']G sel{ (x,x'), (44) 



GadvOr, z) = 26[x', Z(x)]G seU (x, x'). (45) 

The retarded and advanced Green's functions satisfy 
the inhomogeneous wave equation: 



UG{x,x') 



5 D {x~x') 



g 1 / A {x)g 1 / i {x'){d-2)V 



(46) 



where a symmetrization of the delta-density is per- 
formed. Following the above prescription one finds 
the retarded Green's function in terms of the same 
Uk, v (note that the regular term w does not enter): 



Gr, 



1 



2(2tt) 



— e(x'Mx))x 



,^(-i)V^W_, ( _ ff) , (47) 



\m—0 



m 



Here the sum contains terms proportional to deriva- 
tives of the delta function constitute the direct part of 
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the retarded function with the support on the "light 
cone" <t = 0: 



ml 



(48) 



The i>-term is localized inside the light cone and rep- 
resents a tail resulting from the scattering of waves on 
the curvature. 

Using the recurrent relation found for the 
Hadamard function in the previous section, we obtain 
for the retarded Green's function a similar relation 

^k = -2n(d-l)G°+ 2 . (49) 

This relation can be rewritten using the regularized 
delta-functions in terms of differentiation with respect 
to the point-splitting distance £ = e 2 /2, 



S(-a) = lim 8(\<t\-£). 



Obviously we will have 



Wir 



1 dG 



D 



2n(l-d) d£ 



(50) 



(51) 



D. Vector field 

Now consider the massless vector field in D 

dimensions governed by the usual Maxwell action 
quadratic in = V li A v — V„^4 M . The homogeneous 
Maxwell equation 

F" v . „ = (52) 

in terms of A^ in the gauge A^. — reads 

□ A„ - R^A V = 0, (53) 

where R 1 ^ is the Ricci tensor (we use the definition 
as R pv = g p( jR alipv \ Similarly to the scalar case, 
we define the Hadamard function as a solution to the 
equation 

U x G lia {x,x')-R^(x)G 1/a (x,x')=0. (54) 
We look for a solution in the form (D even) 



G 



H 



{2ir) d Vc^- 1 a d - 2 a d ^ 



+ ^1^+^^ + ^ , (55) 
a z a J 

where all the quantities ui (x, x 1 ), v tia {x 1 x 1 ) are two- 
point tensors regular in the coincidence limit. Now 
the lowest order term is 



(56) 



By the same reasoning as in the scalar case we find 
relations between u^ a in different dimensions and es- 
tablish the differentiation rule essentially equivalent 
to (ED. 



III. RADIATION REACTION 

A. Lorentz-Dirac force 

Consider a point particle of mass m interacting with 
massless scalar and vector fields (q , e being the corre- 
sponding couplings). The particle world- line will be 
denoted as x M = z ,1 (s) where s is the interval. The 
action consists of the field part 



S F = -- 



1 



and the particle term 



(V0) : 



1 , 



-gd D x, (57) 



S p = -too / (1 + q$)\J -z 2 dr -el A^dr, (58) 



where <f> and A 11 are taken at the point x = z(t). The 
field equations read 



□ = 



(59) 
(60) 



while the equation of motion of a particle is 

mo(l + #)^ = -eF» v z v - m Q qW v <j> v , (61) 
where the covariant differentiation along the world 



line is denoted by dot, 



d v <j>, and IP" = + 



z^z" is a projector on the subspace orthogonal to the 
world line. The natural parameter t is assumed, so 
that i 2 = —1. 

To get the radiation reaction force we insert at the 
right hand side of the Eq. (|tIT|) the retarded solution 
of the wave equations (|59I60[) 



4>{x) = I G iet (x,x')p(x') y/-g(x')dx', 
A»(x) = J G^^J^^/^gJ^dx', 
with the sources 

p(x ) = m q I ^= dr. 



nx>) = ei^'- z ^ a{T) d T. 



-<] 



(62) 
(63) 

(64) 
(65) 



The resulting Lorentz-Dirac self-force will be by def- 
inition the right hand side of the equation of motion 
in absence of external fields 



mo'z^ = /f- 



LD> 



(66) 



where the mass parameter at the left hand side is con- 
stant (contrary to another frequently used in presence 
of the scalar field definition of the effective mass in 
([5])). The Lorentz-Dirac force consists of a scalar part 



f^ = -m qU^(^ + ^z v ), 
and a vector part 

P = -eF M z v 

Jem oi v" ' 



(67) 



(68) 



6 



According to the above decomposition of the re- 
tarded Green's function into the direct part and the 
tail part, we will have a similar decomposition for the 
Lorentz-Dirac force. Then for the direct part we ob- 
tain the following differentiation rule relating the val- 
ues of the force in twice neighboring dimensions: 



rfi dir 
JD+2 



1 



,£) 



D - 1 



d£ 



(69) 



The limit £ — > +0 has to be taken after the differ- 
entiation. The direct force is due to the light cone 
part of the retarded Green's function. This is not the 
full local contribution to the Lorentz-Dirac force. An 
additional contribution comes from the differentiation 
of the theta function in the tail term v9(—a). In the 
scalar case this contribution vanishes, but in the elec- 
tromagnetic case an extra local term arises: 



•Hoc 



QVua]z u - [v va ]Zu 



(70) 



where the coincidence limit [v va ] depends on the di- 
mension. In the cases for D = 4 it is given by (|92|) . 
The remaining contribution from the tail term will 
have the form of an integral along the past half of the 
particle world line. 



B. Divergences 

The direct part of the Lorentz-Dirac force contains 
divergences. To separate the divergent terms one can 
use the decomposition of the retarded potential sug- 
gested in the case of four dimensions by Detweiler and 
Whiting [25l | . In higher even dimensions we can follow 
essentially the same procedure. We define the "singu- 
lar" part Gs of the retarded Green's function as the 
sum of the symmetric part ( "self" ) and the tail func- 
tion v as follows 



G s (x,x') = G sc u(x,x') + 



1 



^v(x,x') 



4(27r) d " 

= Gself dir(z, x ) + — -p-^ fa x') 0(a). (71) 

Here the direct part of the self function means its part 
without the tail i>-term. Taking into account Ov = 0, 
it is clear that Gs satisfies the same inhomogeneous 
equation as G sc if . The u-term in the second line of the 
Eq. (|7ip is localized outside the light cone. Therefore 
the corresponding field (for instance, scalar), at an 
arbitrary point x will be given by 

Tadv 

(f> S (x) = sc if dir + 4 (2^)d-l J V ( X ' Z ( T )) dT ' ( 72 ) 

Tret 

where the retarded and advanced proper time values 
T re t (a:), T a dv (a;) are the intersection points of the past 
and future light cones centered at x with the world 
line. 

The remaining part of the Green's function 

Gji(x,x') = G iet (x,x') - Gs(x,x') (73) 
satisfies a free wave equation and is regular. 



Let us investigate divergent terms in the Lorentz- 
Dirac force and show that they can be eliminated 
adding counterterms into the action (in flat space this 
was demonstrated for D = 6 by Kosyakov |20j an 
generalized to arbitrary D in 26]). Inserting (l72l 
into the Lorentz-Dirac force defined on the world-line 
x — z(t), we observe that the integral contribution 
from the tail term vanishes and so the divergent part 
is entirely given by <^ S eifdir- For the delta- functions 
we use the point-splitting regularization ([50]) . All di- 
vergent terms will then arise as negative powers of e. 
To prove the existence of the counterterms we con- 
sider the interaction term in the action substituting 
the field as the S-part of the retarded solution to the 
wave equation. In the scalar case we will have: 



1 



S s = - / G s (x,x')p(x)p(x')y/g(z)g(z>) dxdx', (74) 

where a factor one half is introduced to avoid double 
counting when self-interaction is considered. Substi- 
tuting the currents we get the integral over ch~, dr': 

Ss = \m 2 q 2 n J G s (z(T),z(r'))dTdT'. (75) 

Since the Green's function is localized on the light 
cone (by virtue of (|T2"|) ). we can expand the integrand 
in terms of the difference t — t — t' around the point 
z(r): 

SS ~ I dT I^ BH {T)S(k) ^ ~ £2)tldL (?6) 

k.l 

Here the coefficients B^z) depend on the curvature, 
while the delta-functions are flat: S^(t 2 — e 2 ) (the 
derivatives with respect to the argument are under- 
stood). By virtue of parity, the integrals with odd I 
vanish, so only the odd inverse powers of e will be 
present in the expansion. Moreover, if we know the 
divergent terms in some dimension D, we can obtain 
all divergent terms in D + 2 except for 1/e term. The 
linearly divergent term corresponds to I = 2k. The 
integral is equal to 



S (k)(t 2 



; 2 )t 2k dt 



(-l) fc (2fc -1)111 



In four dimensions this term is unique. Applying our 
recurrence chain we obtain the inverse cubic diver- 
gence in six dimensions and calculate again the lin- 
early divergent term. Thus in D = 2d dimensions 
we will get d — 1 divergent terms from which d — 2 
can be obtained by the differentiation of the previous- 
dimensional divergence, and the linearly divergent will 
be new. This linearly divergent self-action term in the 
action will have generically the form 



S. 



(-1) 



E 

k=n-2 



(-l) k (2k- 1)!! 
2^+i 



B k ,2k(T)dT. 



Here the coefficient functions are obtained taking the 
coincidence limits of the two-point tensors involved in 
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the expansion of the Hadamard solution. They actu- 
ally depend on the derivatives of the world line em- 
bedding function z(t) as well as the curvature terms 
taken on the world line: 

Bk,2k(r) = B kt2k (z, z, R(z(t)), R fiv (z(r)), ...)■ 

Finally we have to rewrite the expression in the curved 
reparametrization-invariant form. This amounts to 
the replacement 



z(t) — > Dz = 
z(t) D 2 z, 



1 dz 



\l-z 2 dr 



dr — > \/ —z 2 dr, 



(77) 



where D is the covariant derivative along the world 
line. 

The vector case is technically the same, now one 
has to expand in in powers t in the integrand of 

1 



S S = r / GL (*> WW (x')y/g(z)g(z') dx dx' . 



From this analysis it follows that in any dimension 
the highest divergent term can be absorbed by the 
renormalization of the mass as in the generating four 
dimensional case. To absorb the remaining d — 2 di- 
vergences one has to add to the initial action the sum 
of d — 2 counterterms depending on higher derivatives 
of the particle velocity. 



cone and a tail term: 



G iet = —0[E(x),x'] A^ 2 S(a)-v6(-a) . (78) 

47T L 

The retarded solution for the scalar field reads 



Trot (a:) 



(j>ret(x) =m q 



-A 1 / 2 S(a) +v6(-a)] dr' . (79) 



Differentiating this expression we obtain <fi u on the 
world line: 



4>v{z(t)) =to <7 



-A l ' 2 8'{ 



- A)l 2 5{a) - v5(a)a u + v v ] dr', (80) 

where integration is performed along the past history 
of the particle. All the two-point functions are taken 
on the world-line at the points x — z(t) ("observa- 
tion" point) and z 1 = z(t') ( "emission" point). 

To compute local contributions from the terms pro- 
portional to delta-functions and its derivative, it is 
enough to expand the integrand in terms of the dif- 
ference s = At = t — t' around the point z(t). The 
Taylor (covariant) expansion of the fundamental bis- 
calar <j(z(t), z(t')) is given by (27| : 



o(z(r),z(r')) 



E^^r)(r 

fc=0 



(81) 



C. Finite reaction terms 

To obtain local part of the finite radiation reaction 
terms one can used the recurrent procedure from D 
to D + 2 dimensions. So if one performs expansions 
in terms of t = r — t' in the integrals up to suffi- 
ciently high order, one can find the generating expres- 
sion to get finite reaction terms in higher dimensions. 
In addition one obtains the integral tail term which 
is essentially the same (in terms of function v) in all 
dimensions. 



IV. FOUR DIMENSIONS 

Here we apply our formalism to rederive the ra- 
diation reaction force in four dimensions. This ap- 
proach is much simpler than the non-local DeWitt- 
Brehme type calculations (integration of the stress- 
tensor flux over the world-tube surrounding the par- 
ticle) and gives a more transparent understanding of 
the renormalization involved. Brief account of our ap- 
proach in four dimensions was presented in [19j . Cal- 
culation in four dimensions may also serve as a start- 
ing point for higher-dimensional problems using the 
recurrent relations found above. 



A. Scalar force 

In four dimensions the scalar retarded Green's func- 
tion contains a single direct term localized on the light 



where D is a covariant derivative along the world-line: 
Da = a — <r n z 



D 2 a 



cr Q( gi Q i' 3 + a a z a , 



etc. Such an expansion exists since the difference s = 
t — t' is a two-point scalar itself: this is the integral 
from the scalar function J (— z 2 y/ 2 ds, along the world- 
line from z(t) to z(t'). Taking the limits and using 
z 2 {t) = — 1, we find: 



a(s) = - Y -z 2 (r)-+0( S 5 ). 



(82) 



To obtain an expansion of the derivative of a over 
z^fr) one can expand <7 M (r, r — s) in powers of s. This 
quantity transforms as a vector at z(t) and a scalar 
at z(t'), this 



*•-+ z 



-)+0( S 4 ), (83) 



where the index /i corresponds to the point z(t): o'' = 
da(z, z') /dz^. Recall that the initial Greek indices 
correspond to z(t'). The expansion of S(—a) will read: 



5(-a) = S(s 2 /2) + s 



24 



S'(s z /2) + .. 



(84) 



where the derivative of the delta- function is taken with 
respect to the full argument. Since the most singular 
term is A 1 / 2 5'(a)a v , the maximal order giving the 
non-zero result after the integration is s 3 . (Note, that 
in order to use our dimensional recurrent relations to 
obtain a reaction force in higher dimensions we should 
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perform an expansion up to higher orders in s.) Thus, 
with the required accuracy, S(— a) = 2S(s 2 ), and all 
the integrals for the delta-derivatives are the same as 
in the flat space-time. This allows us to use the same 
regularization of the delta-functions with double roots 
S(s 2 ) by the point-splitting. Expanding the biscalar 
A 1 / 2 and its gradient at z we have: 



A 1 / 2 = 1 + —R„ T z T z°, 



c^A 1 / 2 = -R VT i r . (85) 
6 

Combining all the contributions we obtain finally for 
the field strength: 

( I 1 ... 1 1 

<p v = m Q q \ —z v - - z v --R vr z + —Rz v - 
\ le 6 o 12 



(86) 



where the first term is divergent. Note that the terms 
— 1/6 R\ p z x z p z v + 1/12 Rz v are annihilated by the 
projector IP 1 '. The retarded field 4* itself is also sin- 
gular on the world-line: 



(z) = 1- m q 



vdr'. 



(87) 



Collecting all the orders (|86ti87|) . we obtain the fol- 
lowing expression for the scalar part of the self-force 
(E73: 



2 2 

m q 



( ~ Z u +^R„ T Z T 



v,, dr' 



2e 



vdr' 



The symmetric Green's function is given by 

while the retarded one is 

G Tet (x,x') = 29[T.}G self (x,x'). 
The retarded vector-potential is 

A?(x)=-4*e J G^,z(r'))z a dr' = 

Sret(x) 



(93) 
(94) 



=e J [-u m 6(a)+ Vlia 6(-a)]z a dT. (95) 

— OO 

The field strength on the world-line x — z{s) will be 

s 

— OO 

+v m a u S(a) -{(j, < — >v}\z a dT', (96) 



so the Lorentz-Dirac force will read 

/&,(*) = -eFt(z(s))i v (s). 
Substituting all the expansions 



(97) 



g va (r, t - s)z»(T)z a (T - s) = -1 - z 2 - + 0(s 3 ) 

u v ^z v (r)z a (r -s) = ~R» T * T + 0{s 2 ) 

u\. v z v (r)z a {r -s) = -R vr z v z T z» + 0(s 2 ) 
6 

s 2 

u va a^ l z v {T)z a {T -s) = -sz" + —z"+ 
+ s 3 (~l r -±Rx»^z» - + 0(s 4 ) 



B. Vector force 

Now consider the vector contribution. The 
Hadamard function in four dimensions is 



G 



H p,a 



1 ( u fia , 

1 v„ n In a + w„ 



(89) 



{2tt) z V a 

where u^ a (x, x'), v tia (x,x') and w lia (x,x / ) are bivec- 
tors. One find u MQ = g^A 1 / 2 . We look for power 
expansions 



v« = E <"> n > w ^ = E O n > ( 9 °) 

n=0 n=0 

we find the following equation for 

= -Uu lxa + R ll u u va . (91) 

Substituting the covariant derivative of the two-point 
tensors, we obtain the following coincidence limit for 



1 



l«£«'J = L«M«J = 2 V ( R <*f> ~ ^9apR ) • (9 2 ) 



v? a o. v z v (j)z a (T - s) = -szP + z^s 2 

1 , 



+ s 3 ("^i" - -^i^i" - \ A + 0( s 4 ) 
v va a»z v z a = S -z»R av z a z v + ^Rz" + 0(s 2 ) 
v» a o, v z»i a = - S -R»z a + ^Rz» + 0(s 2 ), (98) 
after integration we find 



-^ + W»(\z v+ \R va z» 



T 

+ i"(T) J {v^-V va ^)z a {T')dT' 



(99) 



C. Renormalization and total force 

Combining the divergent parts of the expressions 
(|99)) and ([87)) into the unique mass-renormalization 
term, we get 



m-m = -^-(e 2 - mlq 2 ). 
2e 



(100) 
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Clearly, this is just the flat-space result. Note that 
under the condition: 



mo\q\ 



(101) 



the model is free from singularities and does not re- 
quire mass-renormalization. The final form of curved 
space Lorentz-Dirac equation in this case will read: 



m(r)^ M = e 2 



W z, 



1 



v,, dr' 



+ Z»(T) (v»„-V^)z a (T>)dT' 



, (102) 



where 



The defining equations are: 



□A 1/2 - 2« /1 o-" - -it | 2 



= 0, 



(106) 



A cr^ 

. \u + 2v^ + v[A ^— ) + aUw = 0, (107) 



Dv = 0. 



(108) 



Note that the representation of the Hadamard func- 
tion in inverse powers of a is non-unique, but this 
amounts only to a redefinition of the regular part w 
which does not contribute to the retarded solution 



4>T{x) 

where 



mq / G"*(x,x')j(x')da/, (109) 



m(r) = m + / vdr' 



is the r-dependent "mass". 

This result coincides with the sum of the DeWitt- 
Brchme-Hobbs and the corresponding scalar equa- 
tions [17| which was previously obtained within the 
DeWitt-Brehme non-local approach. 



V. SIX DIMENSIONS 

The next even-dimensional case is D = 6. The vec- 
tor reaction force in the flat space-time was previously 
given by Kosyakov 20] . Here we give the curved space 
treatment and also add the scalar field. 



A. Scalar force 

Using the general formalism of Sec. [XT] we can 
present the Hadamard Green's function of the scalar 
field in six dimensions as follows 

1 /A 1 / 2 u \ 
G H = J - f (-^ + - + v\no + wy (103) 

The corresponding retarded Green's function will be 
0(£ 



8ir- 



-(A^S'ia) - uS(a) + v0(-of), (104) 



so that the retarded solution is 

Tret (x) 



<t>ret(x)=l J (A^ 2 6'(a)-uS(<j)+v9(-a)) dr', (105) 



where a — a(x, z(t')). 

We have to expand the two-point scalars A 1//2 (x, z), 
u(x,z), v(x,z) around x and then set x = z(t), z — 
z(t'). To get all the local terms explicitly it is sufficient 
to expand u(x, z) up to the second order in keeping 
only the leading term (coincidence limit) for v(x,z). 



8vr 2 



u ;v S(a) - vS(—a)a v + v v 6{—a) 



(110) 



The subsequent calculations consist in expanding two- 
point functions into powers of a, which are now needed 
up to higher orders than in the four-dimensional case 
since the Green's function has higher degree of singu- 
larity. In particular, we get 

A 1 ' 2 =1 + -^firfc/V" - ijS oS; ,,7'V<V+ 



( — 

V288 



60 



(111) 



Also, an expansion of the argument of delta-functions 
will be 



a = la»a„ = I ( -r 2 - ^z 2 r 4 + ifz z)r 5 



2 p 2 V 12 
5{a) = S(-r 2 /2) + S'(~t 2 /2) 



12 



24 



(112) 



Omitting further calculations, we present the final 
result for the scalar Lorentz-Dirac-DeWitt-Brehme 
force: 

m z" = m 2 q 2 (/£ t + div + fin + f^ n ) , 
where the "flat-space" part is given by 



/flat q £ 3 Z + _ 



1 I 1 



16 



-Cz z)z fi - IF" 
9 V ' 



24 " 

-z 2 "z v ^\ 4 5) 

9 45 " 



(113) 



It has the same form as in the Minkowski space, but 
the dots denote covariant derivatives along the world- 
line. It includes two divergent terms proportional to 
e~ 3 and £ -1 , and it is the only part of the total force 
which survives in the flat space limit. 



10 



The divergent part induced by the curvature con- 
tains only e~ l term and is given by 



curv div 



1 



—W v \2R va i a + 2R ua . z a z 

R a /3; v z a iP + R a pz u z a z p + R-U + Rz v 

(114) 



The leading divergence e 3 is absorbed by the renor- 
malization of mass, 



m = mo 



2 9 

mm 

6e 3 ' 



(115) 



it is not affected by the curvature. To eliminate the 
second divergence £~ 1 one has to introduce the fol- 
lowing counterterm into the action: 



K 

72 



3 -(D*zf+ 
(R- R afi Dz a Dz^)\yJ~z^dT, (116) 



where Ko is a new bare coupling constant, and D = 
z 2 -D/dr is the reparameterization-invariant co- 
variant derivative along the world-line. The first term 
in S c is a covariantization of the one obtained in the 
flat space [20] , while the second Ricci-dependent term 
is new. Note that it vanishes in the region of space- 
time free of matter sources. The bare coupling con- 
stant k° entering the counterterm is renormalized to 



(117) 



Non-vanishing k leads to rigid particle dynamics in- 
vestigated in a number of papers (see fTTI j for the ref- 
erences). 

The local part of the finite force induced by the 
curvature depends both on the Ricci and the Riemann 
tensor, it does not vanishes therefore in the vacuum 
region of the space-time: 



— Tpf 

J curv fin 



R P u a RppT~{ RuaRf3j RaRppvi 20 ^"ot'tM + gg Raf3;v~f ) Z Z^ Z 1 



1 



1 



360 



Rj* RotprX n on R ^ R\vpa. a( , R 



180 



60 



— R^R a \ + —R va R + EiR»a i 



a Rol(5 Z %P Zv Rvol Z ~\~ oa R Z v n Ra.fi Z Z^ Z v Roc {3, jZ Z^ Z^ Zy + Rvol'z Z 

R va , z Q i /3 - ^ R va , pz a z? + ^ R aP , u z a z^ 



1 



•a • 8 • -y ■ 



1 



(118) 



In the case of the geodesic motion the last four lines 
in this expression vanish. The non-vanishing part is 
quadratic in curvature. 

The non-local tail term is given by 



" tail 



Uu^ J Vvdr'+z" J vdr'\. (119) 



The retarded potential will be 



TVetO) 



+ v» a 8(-a) z a dr'. (122) 



B. Vector field 



Omitting the details of the calculation we present the 
final result in the form 



In six dimensions the Hadamard function is 



= e 2 (j£ t + /^ urvdiv + jf urvfin + O , (123) 



while the retarded Green's function reads: 



(120) 



G 1 . ct = -—6(x',Z(x))x 
8ir z 

x A 1 / 2 g" a 6'(a) - u" a S(a) + v ua 6(—a)] . (121) 



where again the "flat" part is the covariantization of 
the flat-space expression 



6e 3 



16 



8 " 



+ 2 -(z-zy^-U^(±z^ + h 2 A. (124) 
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The divergent term is 
5 



J ci 



curv div 



72e 



IP" \2R va z a + 2R va ,gz a z 



- R aP , v z a z 15 + R a pz v z a i p + -(R- t „ + Rz v ) 

(125) 

The leading divergence is absorbed by the mass renor- 
malization, while the non-leading one is eliminated in- 
troducing the following countertcrm: 

- (R + 5R afs Dz a Dz^))y/^PdT. (126) 



The finite part of the self- force /^ urvfin can be cast 
into the sum of the following parts: a) the quadratic 
in curvature terms: 



1 



7 



ft* TF" y R R z° ' 7 7 ' 7 P A R R 7 a 4- 

^quad — 11 \ gg-'W-nrp* Z Z f ^ 11 IX V(T Z f 

I 29 

n napar ■ i ___ E> T) aT ^ 



1080"' , T 
1 

540 



4320 



(127) 



b) the Ricci-terms: 



27 



27 



~R ar z' T z T "z +^R ua . T z a z T + -LnR ury z a 
R z T z a —R z a 'z T 



1 



1 



^^Rva\pT^ Z Z^ gg-Rpcr; VT% % ^ J • (1^8) 



c) the Ricci-scalar terms: 



:i „ ... i 

54 ~ " 36" 



1 

72" 



d) the Riemann term 



JRi 



Riem g 1 ovt ** 



r z z z 



(129) 



The appearance of the Riemann tensor in the finite 
local force in a new feature in six dimensions, which is 
absent in the scalar case. This term is non-zero in the 
vacuum space-time, but it still vanishes for geodesic 
motion. 

Finally, we observe that the tail term has the same 
form as in four dimensions. 



VI. CONCLUSIONS 



In this paper we presented a purely local calcula- 
tion of the radiation reaction for a point particle in 
curved space of an even dimension D ^ 4. The pos- 
sibility of such a calculation means that scattering of 
radiation on curvature in the vicinity of a charge is ir- 
relevant for the reaction force. Clearly the scattering 
can substantially modify global properties of radia- 
tion. In particular, emitted waves can be reflected 
back and act on a particle once again, but within our 
approach this should be rather interpreted as an ac- 
tion of the external field, than a modification of the 
reaction force. Note, however, that such a description 
may not be natural if the problem is treated globally in 
terms of the field modes defined on the full manifold. 
In particular, for space-times with reflection of mass- 
less fields from infinity (like anti-de Sitter) the action 
of modes (subject to reflection conditions) on a par- 
ticle will include the effect of scattering, thus leading 
to a global definition of the reaction force. In terms of 
the above analysis this means that there are several 
geodesies connecting two points of space-time. Our 
approach is valid only under the assumption that this 
is not the case. It excludes, therefore, many phys- 
ically interesting situations, which require a special 
global analysis. For an example of the global defini- 
tion of the reaction force see [28[ where the radiation 
reaction problem was considered for particles moving 
in the Kerr space-time, (for the case of a cosmic string 
see [H,[30|). 



Another point worth to be discussed is that of 
elimination of divergences. Classical renormalizabil- 
ity of the equations for point particles interacting 
with massless scalar and vector fields in flat space 
time is well-known for D — 4. In higher dimensional 
Minkowski space it also holds modulo introduction 
of the higher derivative counterterms into the action 
(Til . I20I ]. Here we have demonstrated that this re- 
mains valid in the general curved space-time of an 
arbitrary dimension for a minimally coupled linear 
scalar and vector fields. The curved space countert- 
erms are not simply the covariant generalizations of 
the flat space ones, but generically include curvature- 
dependent terms. An explicit example was given for 
the D — 6 case. One observes that the renormal- 
ized higher-dimensional reaction force contains not 
only higher derivatives of the particle velocity, but 
also higher derivatives of the Ricci tensor, as well as 
the Riemann tensor, a feature which is not present in 
four dimensions. 
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